Quantum structural phase transition in chains of interacting atoms 
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A quasi one-dimensional system of trapped, repulsively interacting atoms (e.g., an ion chain) 
exhibits a structural phase transition from a linear chain to a zigzag structure, tuned by reducing the 
transverse trap potential or increasing the particle density. Since it is a one dimensional transition, 
it takes place at zero temperature and therefore quantum fluctuations dominate. In [Fishman, 
et al., Phys. Rev. B 77, 064111 (2008)] it was shown that the system close to the hnear-zigzag 
instability is described by a 4>'^ model. We propose a mapping of the (j)* field theory to the well known 
Ising chain in a transverse field, which exhibits a quantum critical point. Based on this mapping, 
we estimate the quantum critical point in terms of the system parameters. This estimate gives 
the critical value of the transverse trap frequency for which the quantum phase transition occurs, 
and which has a finite, measurable deviation from the critical point evaluated within the classical 
theory. A measurement is suggested for atomic systems which can probe the critical trap frequency 
at sufficiently low temperatures T. We focus in particular on a trapped ion system, and estimate 
the implied limitations on T and on the interparticle distance. We conclude that the experimental 
observation of the quantum critical behavior is in principle accessible. 
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I. INTRODUCTION 

The structural transition from a linear chain of repul- 
sively interacting particles to a planar configuration, in 
the form of a zigzag structure, has been often discussed 
in theoretical studies on atomic and condensed matter 
systems. Examples include electrons in nanowiresi, ul- 
tracold dipolar gases^, vortex lines in Bose-Einstein con- 
densates^, and ion Coulomb crystals in traps-"^. Specif- 
ically in ion Coulomb crystals this transition has been 
experimentally observed and characterized^i""— , thereby 
determining a phase diagram of the ionic structures as a 
function of the trap aspect ratio and of the mean inter- 
particle distance^. 

Theoretical studies demonstrated that the transition 
from the string to the zigzag chain is associated with a 
symmetry breaking. More specifically, if the Hamilto- 
nian is invariant under rotation about the string axis, 
it is the rotational symmetry around the chain which 
is broken in the zigzag phase. If instead the motion 
of the interacting particles is confined to the plane, the 
zigzag phase breaks the symmetry by refiection about 
the chain axis. The phase transition predicted by the 
classical theory is hence second-orde r 1 ^ ■ 1 ^ . It can be de- 
scribed by a Landau model, of which critical values and 
exponents are well knowni^, and where the soft mode is 
the zigzag mode of the linear chain - namely, the trans- 
verse mode with the shorthest wavelength - which drives 
the instability and determines the new structure. The 
corresponding Ginzburg-Landau equation in the contin- 
uum limit and in presence of damping was reported in 
Ref. [ij, where special focus was given on creation of de- 
fects when quenching the value of the transverse trap fre- 
quency across the mechanical instability. Further works 
discussed the structural transitions when the external po- 
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FIG. 1: (color online) Sketch of the phase diagram exhibiting 
a linear-zigzag transition, according to the mapping to the 
ID model of a quantum Ising transition. Here, T is the tem- 
perature of the sample and the dimensionless parameter e is 
tuned by the confining potential or the interparticle distance. 
The value e = corresponds to the value of the transverse fre- 
quency, where the linear chain is classically unstablelS. The 
quantum critical point, at > 0, separates the linear chain 
from the zigzag phase at T = 0. For < e < Ec quantum 
fiuctuations dominate, and the crystal is in the linear (dis- 
ordered) phase. The dashed lines indicate the boundaries of 
the quantum critical region, where thermal fiuctuations dom- 
inate. 



tential is not harmonioii., showing that the statistical me- 
chanics of the system at the instability may be profoundly 
modified. The results of these studies are strictly valid in 
the classical regime, when the thermal fluctuations can 
be neglected. 

In Ref. [1^ we addressed the question, whether quan- 
tum fluctuations significantly modify the location of the 
critical point and the behavior in its vicinity, focussing on 
ion Coulomb crystals. This question draws on numerical 
studies in low dimensional dipolar systems, which showed 
a significant contribution of quantum fluctuations at the 
transition from a string to a zigzag order^^. In this article 
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we provide the detailed calculations at the basis of the 
results presented in Ref. [l^ and extend the treatment 
to interacting atomic chains, with repulsive interaction 
potential scaling with with a > 1. We focus on the 
case where the motion of the particles is confined to a 
plane. As argued below, the string-to-zigzag transition 
is indeed a quantum phase transition analogous to the 
ferromagnetic transition of an Ising chain in a transverse 
field^^. A phase diagram summarizing the behavior in its 
vicinity is depicted in Fig. [TJ Note that a similar map- 
ping to the Ising model was argued in the context of elec- 
tronic systems in quantum wiresi. However, the validity 
of the assumed equivalence between a continuous field 
theoretical model and a discrete {Z-i) field-theory cannot 
be justified on the basis of classification into universality 
classes below three space dimensions^*. 

In the present study we develop a quantum field the- 
oretical approach, which allows a direct mapping to the 
quantum one-dimensional Ising model under plausible 
assumptions. On this basis we determine the quantum 
critical point at which the transition occurs and iden- 
tify the experimental parameters, for which it can be 
distinguished from the value predicted by the classical 
theoryJ^. In particular, we relate the critical value of 
the transverse trap frequency vt to the ratio of the typ- 
ical kinetic energy scale \Jk and interaction energy Up. 
Defining a parameter e — 1 — /v^, which describes 
the deviation of the transverse confinement from its sta- 
bility point, we find that the critical point is given by 
Ec ~ {Uk /Up)"^^^ . A measurement is suggested, which 
can demonstrate transverse quantum fluctuations in a 
long chain of trapped ions (composed of hundreds of ions) 
and probe the different regimes in the phase diagram. 

This paper is organized as follows. In Sec. |n]we de- 
rive a quantum field theory for the zigzag phonon mode. 
The mapping to the quantum Ising model is described 
in Sec. IIIII In Sec. IIVI we derive expressions for the 
Ising model parameters, and in particular their critical 
values, in terms of microscopic parameters of the inter- 
acting atoms system. In Sec. |V] we discuss an experi- 
mental realization in trapped ions systems, and specify 
the conditions under which the quantum phase transi- 
tion is likely to be observable. Finally, our conclusions 
and outlook are summarized in Sec. I VII 



II. QUANTUM FIELD THEORY FOR THE 
ZIGZAG MODE 

We start by considering a linear array of particles with 
mass m, which are confined in the x — y plane. The 
system is at low temperature, and the particles at equi- 
librium are aligned along the x-axis and equidistant, with 
interparticle distance a such that their equilibrium posi- 
tion is xf^ = ja and yj^^ = 0, while the motion along 
the z-axis is assumed to be frozen out. Periodic bound- 
ary conditions are assumed at the chain edges and one 
ion of the chain is assumed to be pinned. The particles 



can oscillate around the equilibrium positions in the x 
and y direction, where they are confined by a harmonic 
potential of the form 

1 ^ 

y, = i™.,^^y|; (1) 

here yj is the displacement of the j'th atom from the 
chain axis and vt the trap frequency. Finally, they inter- 
act by means of the repulsive power law potential, 

^"*"^£[(^*-^^-)^ + (y.-%)^]"/^ (2) 

where {r^ = (xj,yj)} denotes the atomic position in the 
plane, is the strength of the interaction and a> 1. 

This model includes the Coulomb interaction when 
a = 1 and the dipolar interaction when a = 3, assuming 
that the dipoles are polarized along the z axis, such that 
in the x — y plane they experience an isotropic repulsive 
force. We also note that in the case of the Coulomb in- 
teraction the strength is such that A = Q, where Q is the 
charge of the particle. In this case the array can result 
from the self-organization of the atoms, since quasi-long- 
range order can be assumed as long as one particle in the 
chain is pinned. The configuration, in which the particles 
are aligned along the x-axis with uniform interparticle 
distance, could be realized in a ring trap at ultralow tem- 
peratures under the condition that one ion is pinned^, or 
in a anharmonic linear Paul trapi^. It also corresponds 
with good approximation to the ions distribution at the 
center of a long chain in a linear Paul trap^^i^i^Si^i. 

For a > 1, on the other hand, particle-particle correla- 
tions decay at T = with a power law dependence on the 
distance^^, and long-range order may be only assumed 
in presence of an additional periodic potential along x, 
which localizes the particles at the positions For 
dipolar atoms or molecules with an electric dipole mo- 
ment D, this configuration could be realized by means of 
an optical lattice^^ with periodicity a. Here, the dipole- 
dipole interaction can be made repulsive in the x — y 
plane when a static electric field along the z axis aligns 
the dipole. Then, the strength of the interaction is such 
that A(x D. 

We assume temperatures such that the atoms crys- 
tallize and their motion is well approximated by har- 
monic vibrations about the equilibrium positions when 
the configuration is mechanically stable. While the linear 
chain is stable against longitudinal displacements from 
the equilibrium positions, a soft mode in the transverse 
direction drives a structural phase transition of the lin- 
ear chain into a zigzag configuration, as shown in Ref. [T2I . 
This mode has wavelength A — 2a, such that in the clas- 
sical harmonic crystal the corresponding displacements 
are described by the function yj = {—lyyo with yo the 
amplitude of the oscillations. We therefore focus on the 
dynamics of the transverse phonon modes close to the 
instability and expand Vint to fourth order in \yj\ <C a. 
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Moreover, we introduce the field which describes a 
transverse excitation close to the instability and is slowly 
varying with respect to the length scale a, i.e. 

0. = (3) 



Using a gradient expansion to leading order—, we obtain 
the effective potential 



Terms describing the coupling between the axial and 
transverse modes have not been reported, because they 
give rise to higher order corrections in the gradient ex- 
pansion. The latter is defined by the small parameter 
Ska ^ 1, with dk — \k — ko\ the typical deviation be- 
tween the wave vector k of the excited mode and the wave 
vector fco = 7r/a of the soft mode. An extensive discus- 
sion on the derivation of the effective potential for the 
soft mode can be found in Ref. [H. The ordered, zigzag 
phase is established when acquires a finite, constant 
expectation value, as we show below. To leading order 
in a gradient expansion, the potential in Eq. (U) can be 
mapped to the dynamics described by a potential of the 
form 



For e > 0, the local potential Vq [Eq. (jS])] has minima at 



1 A 



N 



, (5) 



where the interaction is now nearest neighbors. Here Vq 
is a local potential 

Voicj)) = -\m{vl - vl)a^4>^ + ^gaU' , (6) 
and the parameters K, and g are given by 

Ci[a]aA^ 2 4C2[a + 2]aA2 
K = : , — 



a" ma 
8C2[a + 4]a(a + 2)^2 



Q-l-2 



(7) 



with 



CM = 



^ ( \ 



J' 



' E ^ 



A classical theory, which neglects quantum fluctuations 
(i.e., valid in the limit h = Q), predicts a transition from 
a linear to a zigzag chain at T = when the transverse 
confining frequency Vf is reduced below the critical value 
For later convenience we define a dimensionless 
parameter describing the deviation of the transverse trap 
frequency from this classical transition point: 



(8) 



±0 



where the frequency lu is given by the relation 




(9) 



(10) 



and is associated with the curvature of Vq at the top 
of the barrier between the minima of the double well 
potential, as one can observe from Eq. ([5]). 

We note that the mapping to Eq. ([5]) is valid for a 
power-law interaction with a > 1 (for which the rele- 
vant sums Ci[a\ in Eq. ([7]) converge), hence including 
the case of the Coulomb interaction at a = 1. Therefore, 
for the linear-zigzag instability the Coulomb interaction 
effectively belongs to the class of short-range potentials. 
Indeed, one finds that the spectrum of the transverse ex- 
citations close to the instability (approaching from the 
side where the chain is stable classically) behaves as 

a;2 - + {K/m)k^ 

for a > 1 (sec for instance Ref. [l^ for the case a — 
1). This behavior is markedly different from that of the 
long wavelength, axial excitations: for a > \ the axial 
frequency cj ^ ck, with c sound velocity, while for a = 1 
there exists no sound velocity, and w ~ ky/\ log fc p°i^^'^^ . 

These considerations on the structure, which are based 
on identifying the minima of the potential energy, do not 
account for the contribution of the kinetic energy. Ther- 
mal effects are expected to modify the behaviour at the 
critical point. Moreover, even at temperatures T — >■ 0, 
when thermal fluctuations are small, quantum fluctua- 
tions will become relevant. Sufficiently close to the value 
Vc, i.e. for sufficiently small values of the parameter e 
(where the classical zigzag ordering is expected), quan- 
tum fluctuations will induce tunneling between degener- 
ate minima of the potential and are expected to destroy 
the zigzag ordering. A true phase transition will there- 
fore occur at a smaller value of vt, corresponding to a 
quantum critical point, such that vt < Vc] i-e., at Cc > 0. 

In order to investigate the quantum critical behavior, 
one has to introduce a 1-1-1 dimensional quantum field- 
theory for the system. To this end, we write the partition 
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function as 



V(j)e 



-S[4>\/h 



with the Euclidean action 



S[4>] = / dTj2 



ma'^{dr(j)j)'^ 



2'-\r] 4'j+iy 



(11) 



(12) 



zigzag phonon mode in terms of a real continuous scalar 



field 



< 0j (t) < oo. We now express it in a form 



which will allow its mapping onto an effective model for 
the discrete field crj^r) = Sgn[(/)j(T)]. As a starting point, 
we perform the standard procedure of dividing the imagi- 
nary time (t) axis into discrete time steps {rm} separated 
by an infinitesimal interval of size St, with m = 0, . . . , Af 
and M = H/S/St. The partition function is then cast in 
the form 



where /? = l/ksT and /3 -> oo for T — 0. Below we 
demonstrate the mapping of this model to the Ising chain 
in a transverse field, and derive approximate expressions 
for its parameters in terms of microscopic parameters of 
the interacting particles system. This will facilitate the 
study of the conditions under which the phase transition 
can be observable in such systems. 



Z = Tt 



(13) 



III. MAPPING TO THE QUANTUM ISING 
CHAIN 

The low-energy model derived in the previous section 
[Eqs. ([TT|) . describes the quantum dynamics of the 



where T is a transfer matrix describing the time evolution 
from Tjn to T,n+i- Using the notation (j>j — (j)j{Tm), 4>'j — 
4>j{Tm+i), the matrix elements of T are given by 



N 



St 



N 



ma 



St 



^ + V^o(0j) + V^{cj,'^) + \k ((0, - 0,+i)2 + (0^. - 0;.+i)2) 



.(14) 



(15) 



Here the propagator G[(j)j,(j)',\ST) is given by 



tion) therefore acquires the form 



G(0„0';<5t) = 



(16) 



in which the local Hamiltonian Hj — p]/ {2m) + Vq de- 
scribes the quantum dynamics of a particle in the double- 
well potential Vq (</)), where pj is the momentum con- 
jugate to The symmetry of the potential implies 
that the exact eigenstates of Hj consist of pairs of states 
with well-defined symmetry under (p ^ —(p: symmet- 
ric (l^n)) and antisymmetric (|A„)). The corresponding 
energy eigenvalues can be expressed as /i(w„ — Acj„/2), 
h{ujn -\- lS.LOn/2) respectively, where the splitting energy 
h/S.uJn is associated with the overlap of wave functions 
centered in either of the two wells. The propagator 
G{4>, 4>' ] St) (the index j is dropped to simplify the nota- 



G{<^,(^'-St) ^Y.^-'^'"'' [{P'\Sn){Sn\ 



{P'\An){A,, 



r/2| 



(17) 



We now change basis into wave -functions centered at 
the right and left potential minima 0± [Eq. ©]: 



|i?„) ^ -i=(|5„) + |A„)) , |L„)^-i=(|5„ 



\An)) 



(18) 
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Equation p7)) becomes 

G(0,0';<5r)=^e-'^"^ 



(19) 



X ^[{cj)'\Rn}{RrM) + ((/.'|X„)(L„|0)] cosh 
+ [(0'|i?„)(L„|(/)) + (0'|L„)(i?„|0)] sinh 
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We next define the wave functions '^,ix{<f>) = (0|Cti)i 
where C = +, — denotes the isospin R, L, respectively. 
Since these wave functions are not known exactly, we 
implement a variational approach and assume the trial 
function 



'i'nxW = /«(</>) exp 



211 



(20) 



where fn{4') for ?^ > 1 is an oscillatory, symmetric func- 
tion of (j) with a number of nodes increasing with n; 
for the lowest energy states n = we assume fo[(f)) — 
l/i^iy/'^, so that I = lo is a single variational pa- 
rameter. The dependence of 5'n.^(0) on the magnitude 
and sign of (j) can be made explicit using the substitution 
(j) = |(/'|o', which yields 



"^n.xW = /„(|0|)exp 
X {6ac + exp 



(101 



-210100 



(21) 



The first term in Eq. (pij) dominates as long as (/>o > /„. 

To complete the derivation of an effective field the- 
ory in terms of the discrete fields (Jj{T), one needs to 
perform the path integral [Eq. (fTTj) ] over the magnitude- 
field |(/)j(r)|. Due to the Gaussian factor in Eq. (PT|) . the 
integration over = 4>j{T) (for each n) can be written 
schematically as 



; exp 



2/2 



exp{F{(Pj)ST} (22) 



where F{(j)j) encodes the remaining (/>j- dependence, in- 
cluding in particular the interparticle coupling terms 



{K/2){(j)j - (l)j±if [the last term in Eq. If we 

now consider the case where (j)o ^ /„ (which is valid 
deep in the zigzag phase), a saddle-point approximation 
of the integral yields (up to a multiplicative constant fac- 
tor) e^^'^°'^^^^^ . This implies that the partition function 
[Eq. (fT3|) ] can be recast as 



Z = Tr 



(23) 



with Tr = Teff a 2^ x 2^-matrix related to 

the original transfer matrix T by 

TefAWA, W,}] = nicl^oa,}, {00^ ■}] . (24) 
Recalling Eq. (fT4|) . we find that Teff[{(Jj}, Wj}] can be 
expressed as 



Teff[{(Tj},{crj}] 



(25) 



N 



where T/oc[o', cr'] = const x G{(/)oa,4)Qa';ST), and the 
propagator G{(l)ocr, (f)oa' ; St) can be obtained by insert- 
ing Eq. (EI]) in Eq. 1^. 

The final stage in the mapping to the quantum Ising 
chain involves a derivation of an explicit expression for 
the 2 X 2-matrix Tioc, which encodes the single parti- 
cle dynamics. By inspection of Eqs. ((2T|) and (fT9|) it 



is evident that Tioc obeys the symmetry Tioc [c, a'] 
Tioc[~cr, — cr']. It can therefore be written in the form 



const X {Aoa° + Axo'^) 



(26) 



where cr" are Pauli matrices in the basis where \a = ±1} 
denote the eigenvalues of cr^, and is the 2x2 unit 
matrix. This is the crucial point where we formulate the 
symmetry of the problem, which stems from the sym- 
metry relations leading to Eqs. p^ - (IT51) . The general 
expressions for and Ax are 



J 



^0 = E' 



f2(0o) (^coshj 



r smh -i — ^ — 



r 



Ax = ^e--"^V^(0o)(sinh|^^^|+r2cosh 



(27) 



where 



(28) 



Equation (P7|) implies that in general Aq and Ax are 
non-trivial functions of the time interval 5t. Restric- 
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tions on the parameters must therefore be set in order 
to estabhsh a well-defined effective model for arbitrarily 
small St. A simplification can be achieved in view of the 
fact that the sum over n is rapidly converging: a semi- 
classical calculation indicates that a;„ ~ n'^^'^ for large 
n (see Appendix for details), i.e. the level spacing in- 
creases with n. In the regime where the tunnel splitting 
obeys Aw„ <^ w„ (realized far enough from the quantum 
transition so that the wells of Vq are well separated), 
the sum is dominated by the first term n = up to ex- 
ponentially small corrections. Even closer to the phase 



transition, where the tunnel splitting becomes compa- 
rable to the barrier between the wells, the sum in Eq. 
([27l) is dominated by the first few terms. As we show in 
the next section, provided the quartic term in Vo (i.e., 
the parameter 17) is sufhciently large, all the energy lev- 
els h{ujn ± Awn) are given (up to a numerical factor) by 
the same energy scale {h'^g/m'^y^^. Since the differences 
(Aw„ — Awn-i) ~ n~^/^ are typically smaller than Aujn, 
the dependence on n of all terms in the brackets in Eq. 
(P7)) may be neglected. One can therefore approximate 
AuJn for all n by Aw ~ {Hg/m'^y^^. This yields 



J 



Tie 



const X 



, f AluSt 

cosh < 

I 2 



sinh 



AloSt 



l^sinh 



AujSt 



cosh 



AujSt 



(29) 



r 



Assuming in addition that the overlap factor F [Eq. (l28l 
obeys 



r < tanh{Aa;(5T/2} < 1 , 
Eq. dSll) reduces to 



(30) 



Tioc[cr, cr'] 



(31) 



^ , , . AujSt 1 . , f AujSt . J. 
To cosh <^ — - — )■ baa' + smh \ ) d^.-a' 



where Tq is a constant prefactor. In terms of Pauli's 
matrices this becomes 



Tioc ~ To cosh 



AujSt 



(t" + snih i — - — ^ a' 



= Toexp{^}. (32) 

As a result, the expression for Z [Eqs. (|23|) and ([25l) ] 
effectively reduces to 



(33) 



in which Zq is a constant and Sj is the action of an Ising 
chain in a transverse fiel d^^'^^ subject to the Hamiltonian 



TV 



Hi^-J2iJa]cT]^, + ha]) 
i=i 



(34) 



Here, the fictitious exchange coupling J and transverse 
field h are given by 



HAiv 



(35) 



where we have used Eq. ^ for (j>Q. The model described 
by Hamiltonian (1341) is known to exhibit a quantum phase 



transition at h/ J — 1 and T — Q, separating an ordered 
phase at J > h (in our case, the stable zigzag configura- 
tion) from a disordered phase ai J < h (the linear chain) . 
In both phases the spectrum of excitations is character- 
ized by a gap with energy A = 2\J —h\. The correspond- 
ing phase diagram, showing the critical behavior as a 
function of e (and hence vt) and T is sketched in Fig. [1] 
The critical regime {ksT > A) is characterized by uni- 
versal power-law T-dependence of correlation functions 
with the critical exponents of the quantum transitioiiii. 
In the next section we evaluate the gap in terms of the 
atom chain system parameters, and in particular relate it 
to the parameter e [Eq. ([5])]. This will enable an estimate 
of the critical value ec satisfying A(ec) — 0. 

It should be noted that the condition (PH)) was required 
to obtain the form (|32|l for the local transfer matrix, and 
thus the final stage in the derivation of Eqs. (|33 l) . (|34l) . 
Namely, the mapping to the Ising model Eq. (|34|) is valid 
only as long as the lowest energy levels are below the 
barrier in the potential Vq. In this regime, the tunneling 
between the wells is small and hence necessarily h <^ J . 
This is satisfied deep in the ordered (zigzag) phase, in 
which case all our approximations are justified. It does 
not hold near the critical point, which we estimate in the 
next section. However, we note that conformal symme- 
try strongly restricts the number of universality classes 
in f + f-D Moreover, it was argued that the number 
of relevant operators in two dimensions is the same as in 
4 — e dimension^^, therefore the resulting 0* field theory 
describing our model is expected to be in the same uni- 
versality class as the two dimensional Ising model. Based 
on the symmetry properties of the present system (as 
manifested by Eq. (PSI) and the following derivation) it 
is reasonable to assume that a quantum Ising model is 
the appropriate field theory, which correctly describes its 
critical behavior at h ^ J. This hypothesis is further sup- 
ported by numerical studies, e.g., by Barma and Fisher— 
and a later work by Kim, Lin and Rieger— . Under these 



7 



assumptions we estimate the quantum critical point in 
the regime of parameters where classically zigzag order 
is suppressed by quantum tunneling, and find that it be- 
longs to a universality class^'''^'' that differs from the one 
of the classical Landau model. 



IV. ESTIMATE OF THE QUANTUM CRITICAL 
PARAMETERS 



As we show below, the effective coupling parameter 
h/J oi the Ising model [Eq. ([55]) ] can be expressed in 
terms of two dimensionless parameters characterizing the 
system: the first is the ratio of energy scales Uk/Up, 
where 



Uk 



(36) 



is a typical kinetic energy scale of atoms in the chain, 
and 



a" 



(37) 



is an energy scale associated with the interaction. In 
atomic systems, this ratio can be tuned most effectively 
by controlling the density of atoms 1/a, e.g., by means 
of the lateral confinement. The second dimensionless pa- 
rameter e [Eq. ([S])] can in principle be tuned indepen- 
dently by controlling the transverse confining potential 
frequency vt- 

The expression for the fictitious exchange coupling J 
in terms of the above definitions can be obtained directly 
from Eq. Employing Eqs. ©, ^ and §7}, we 

rewrite the parameters as 



K — C\[a]aUp , muj 



2 _ 4C2[a + 2]aC/pe 



8C2[a + 4]a{a + 2)Up 



This yields 



J = CjUp e where Cj 



Ci[a]C2[a + 2]a 
2C2[a + 4\{a + 2) 



(38) 



(39) 



Note that since we focus on the regime e > 0, and by 
definition e < 1, the energy scale J can be tuned contin- 
uously from to ^ J7p. A quantum phase transition is 
expected to occur at ec, dictated by the ratio Uk/Up at 
which J is equal to the "transverse field" energy scale h. 

To get a concrete estimate of h in terms of Uk and 
Up, we perform a variational calculation of the splitting 
energy hAtj. To this end, we consider the trial wave 
functions 



^'±(0) 
X ( exp 



1 



(0 - 0o) 



(40) 



2P 



± exp 



2/2 



for the lowest energy states [i.e., the symmetric and an- 
tisymmetric combinations of Eq. pO|) with n = 0] , and 
search for a minimum of the energies 



E± = 



{±\Ho\±) 
(±\±) 



(41) 



with respect to the variational parameter I. Here Hq is 
the Hamiltonian of a particle of mass m subject to the 
double-well potential Vb, and \E'±((/)) = (0|±). It is con- 
venient to eliminate the parameters of the quadratic part 
of Hq by introducing the length scale A = ^Jh/muj and 
the energy scale huj/2. Using the normalized coordinate 
(j) = {a/X)(j) we write the Hamiltonian as 



- _ d2 



(42) 



depending on a single dimensionless parameter 

The nature of the eigenstates and eigenvalues of cru- 
cially depends on B being small or large. 

The exact minimum condition for E^[l] yields a cum- 
bersome, transcendental equation for the variational pa- 
rameter /. In particular, it includes the exponential fac- 
tors e-'^o/'^, with 00 = (a/A)0o, I = {a/X)l. These 
factors are associated with the overlap of the right and 
left-centered wave functions. However, since we are in- 
terested in a regime of parameters close to the critical 
point (see discussion between Eqs. (|46l) and (|T7)) ). we 
focus on the approximation of a "quasi-flat" potential, 
corresponding to i? 1 and ^q/P ^ 1. Essentially, this 
regime amounts to a dominance of the quartic term in the 
potential, while the barrier separating the two wells is rel- 
atively small. In fact, the parameter B can be written 
a,s B — hbj /{2A.V), where /SV is the difference between 
the energy of the minima of Vb and the top of the bump 
separating them. Therefore, B is the ratio between the 
zero-point energy fiw/2 and the potential energy of the 
barrier. As becomes clear from the following discussion, 
the regime i3 > 1 is relevant for the critical point. In this 
regime also the two Gaussian functions of Eq. pO|) are 
expected to mix, therefore ^g/P < 1. To find an analytic 
expression we take these inequalities to the extreme. We 
find from Eq. (HI 0o — '^'^ leading order in 1/-B, this 
yields a minimum of E-\- [I] at 



Ir. 



A; 



i-mm 7 ''mm 



1 

ZB 



1/3 



(44) 



Both [I] and [I] are computed using Eq. (|^T|) within 
the above approximations, with the value of /^^^ which 
was found minimizing i?+ [I] . We checked that would we 
extremize E- [I] , a value with the same dependence on 
parameters that differs only by 3% would be obtained. 

The resulting energy eigenvalues E± = E±[lmin] are 
separated by a splitting energy Ai? = E^ \l] — E^ [Z] where 

h^^^\h.m"' ■ (45) 
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This sets the scale of the transverse field. 

Recalling Eq. (^5]) for B, we note that the dependence 
on uj (and hence on e) cancels out. We thus get 



1 /3ft^ 

2 V2to2 



1/3 



^ _ /3C2[a + 4]a(a + 2 



1/3 



(46) 



where we have used Eq. (1381) to relate g to Up. It should 
be pointed out that (up to a numerical factor close to 
unity) the above estimate for h coincides with the energy 
scale dictating the higher energy eigenvalues of Hq in 
the semiclassical approximation (formally valid for n 3> 
1): En ~ n'^^'^h [see Appendix, Eq. (|A5|) ]. We hence 
conclude that in this regime of parameters, all energy 
levels are determined by the same energy scale. 

We finally employ Eqs. (jJH) and to evaluate the 
critical value of e for a given Uk/Up: the condition 



dictates 



Up 



2/3 



J 



12(a + 2)4^ 



(47) 

where we have used the approximation Ci [a] ^ C2 [a] ^ 1 
for arbitrary a. For a = 1, a = 3 the numerical pref actor 
is Cc ^ 10. This implies that for given Uk/Up <^ 1, 
a quantum critical point is typically expected at Ec 3> 
Uk/Up. 



V. EXPERIMENTAL REALIZATION 
A. Ion Coulomb crystals 

An interesting system to study the above predicted 
quantum phase transition would be an ion crystal in a 
quasi one-dimensional trap^ '^°i^^'^°i'^^ . For singly-ionized 
alkali-earth metal atoms, the generalized charge A is the 
electron charge {A = e), and the mass of the ions can be 
written as m = UAmp where is the proton mass and 
UA the atomic number. In these systems it is possible to 
tune through the critical point by either controlling the 
transverse confinement frequency ft (which is typically 
in the MHz regime), or the spacing between neighboring 
ions a which is conveniently measured in units of microm- 
eters: a = flo X 1 fim. 

To be able to distinguish between the quantum disor- 
dered phase and the ordered (zigzag) phase, one must 
first make sure that the frequency difference Siy — Vc — i^t 
implied by our estimated Cc [via Eq. ([8])] is not limited by 
the experimentally accessible resolution. To this end, it 
is useful to write in terms of ua and ag defined above. 
Inserting Eqs. (|36p . ([57)) and the numerical values of h, 



e and irip in Eq. (|47l) . we obtain 



(uAao)'^/^ 



(48) 



From Eq. ([8]) , we therefore obtain an upper bound to the 
frequency resolution 



Siy w 10-* 



1 



2(nAao)2/3 



(49) 



which ranges from 100 Hz (for protons with oq ^ 1) to 
1 Hz for heavier ions. This bound has to be compared 
with the time scale of heating in ion traps. In order to 
be able to measure the quantum phase transition, the 
heating and decoherence time scale of the trap, Tdccoh 
should be such that SvTdecoh 1, which leads to de- 
manding conditions for the existi ng trap ping setups, see 



LgtrapE 



for instance the discussion in Ref. 

Another requirement on the experimental setup is the 
possibility to reduce the temperature T below the energy 
scale characteristic of the gap A ~ /i. Using Eq. (PS|) . 
the estimated restriction on T (defined for convenience 
in units of milli-Kelvin) is therefore given by 



T[mK] < 



103C;,(C/i[/p)i/3 



(50) 



To derive a numerical estimate, we again employ Eqs. 
pop. p7p and rewrite the mass and inter-ion spacing in 
terms of ua, oo- This yields 



r[mK] < 0.25 



1 



1/3 



(51) 



implying an upper bound of order ^ O.lmK for protons to 
fiK for Magnesium ions. These values are challenging for 
large crystals, but could be accessed when, for instance, 
only the transverse motion is cooled (Methods for cool- 
ing ion Coulomb crystals to ultralow temperatures are 
discussed in Ref. 33,34). 

We note that the condition on dv can be relaxed in 
the presence of screening with permittivity e > Eq, with 
£0 the vacuum permittivity: in this case the potential 
energy Up decreases as Sq/s, and the upper bound on 
Sv increases with e^/^^ jt is interesting to note that, 
correspondingly, the condition of the temperature be- 
comes more restrictive, as the upper bound decreases 
with e~^/3_ This situation could possibly be realized, 
for instance, in the case of an ion chain embedded in a 
crystal, as in the experiment reported in Ref. [21I . 



B. Dipolar gases 

Another physical system, where the quantum phase 
transition can be experimentally observed, are ultracold 
dipolar gases, such as polar molecule^. Here, the inter- 
action is repulsive provided the dipoles are aligned per- 
pendicularly to the plane where their motion is studied. 
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FIG. 2: (color online) Structure form factor S{ky)/N as a 
function of ky (in units of the transverse length <j)oa/n): (a) 
For an ion string (e < tc). The red dashed curve corresponds 
to the classically stable regime (i), where S{ky) is given by 
Siin [Eq. (I53|) ]: the blue solid curve corresponds to the ion 
string in the quantum disordered phase regime (ii) (Sdis in 
Eq. (I56p with (j)o/l — 1.5). (b) For a zigzag structure, regime 
(iii) (with 00 = 10/). The dashed and solid lines correspond 
to Szz{ky) [Eq. (|55p ] when kx = vrn/a with n even and odd, 
respectively. 



for instance, in a very steep trap which freezes out the 
motion in the perpendicular direction. In Ref. fH a zigzag 
ordering across tubes, containing gases of ultracold polar 
molecules, has been predicted when the distance between 
the tubes was below a certain value. Long-range order is 
found in the presence of an additional external periodic 
potential which localizes the dipoles in an equidistant ar- 
ray in one dimension and could be realized by means of 
an optical lattice. The regime here studied would then 
correspond to the realization of a Mott-insulator state 
of the polar molecules in a one dimensional optical lat- 
tice, where the transverse direction must be confined by a 
dipole trap, whose steepness is changed in order to reach 
several points of the phase diagraro^^. An estimate of the 
parameters, which allow to access the low dimensional, 
quasi-crystalline region, has been presented in Ref. [l6i . 



C. Measurement: structure form factor 

In atomic and molecular systems, the measurement of 
the disordered phase, where quantum fluctuations dom- 
inate in the transverse direction, can be performed by 
means of the structure form factor, here defined by the 
expression 



S{k) 



1 

N 



dx / dx e 



/ ik-(x-x') 



(n(x)n(x')) (52) 



where n{'x.) = n{x, y) is the density of atoms in the 
two dimensional plane. For wave vectors in the limit 
kx ^ kg, where kg — ir/a is the wave number of the 
zigzag, S'(k) — Sf + SoS{kx — fco) in which Sq is propor- 
tional to the squared order parameter, cx (/)q, while 
Sf is the contribution of the fluctuations and is pro- 
portional to the isothermal susceptibility x oi the cor- 
responding Ising system^ Approaching the quantum 
critical point for T ^ A (i.e., for e — > Ec below the dashed 
lines in Fig. 1), x ^ kc— el"'' and 0q ~ |ec— ep'^ in the or- 
dered, zigzag phase (e > Ec), where f3 — 1/8 and 7 = 7/4 
are the exponents of the classical two dimensional Ising 
model. In the critical region T ^ A (shaded area in Fig. 
1), these behave as x ^ ^2 ^ yi/4 respectively, 

where we have used the critical scaling of the correlation 
length ^ - 1/A - |ec - e]'" with v = l21. Note that 
the predicted critical behavior is strictly valid for an in- 
finite system and assumes a uniform density of particles. 
However, in practice it is expected to be approximately 
valid, e.g. in the center of a linear Paul trap of ions, when 
the variation of the density in a large region is negligible 
within a correlation length. 

The distinct regions in the phase diagram (Fig. 1) 
can be identified experimentally also away from criticality 
(i.e., for k that is sufficiently different from kg) with the 
help of light scattering at the appropriate wave vector. 
The structure factor will take a different form in each of 
three cases: (i) the linear (string) chain (e < 0), (ii) the 
quantum disordered phase (0 < e < Ec), and (iii) the 
zigzag configuration (e > €c) for iV ^ 1 atoms. 

In regime (i), where the ion string is assumed, the 
atomic density can be written as n{x,y) — nx{x)ny{y) 
such that Uxix) is the sum of well spatially-separated 
Gaussians, centered in ja and with width Ix <?C a, while 
riyiy) is a Gaussian with width ly — la such that ^ ^ 1. 
The structure form factor is then given by 



1 S{kx J ky) 



with 



5(fc„fc,) = e-'='''/2g-('^'''^)'''/2 



(53) 



(54) 



and corresponds to (5-peaks along the a;-axis at the loca- 
tion of the vectors of the reciprocal lattice, G„ = 27rn/a, 
where the height of the peaks is embedded in a gaussian 
function of width \/2/lx, while the dependence on ky has 
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the form of a gaussian function of width \/2/{la) (see 
Fig. Ha)). 

We next consider the zigzag ordered phase regime (in). 
Denoting by yo — (jioa the transverse displacement from 
the a;-axis, for 0o 3> / we obtain 



kya(j)Q 



(55) 

which provides an interference pattern in the y-direction, 
with a maximum or a zero at fcj, = depending on 
whether the corresponding x component is even or odd, 
as displayed in Fig.[2jb). We note that the appearance of 
an interference pattern for n odd (i.e., the doubling of the 
unit cell) is a primary signature of this ordered phase. In 
both cases, the pattern is characterized by peaks sepa- 
rated by tt/ a0o . These expressions are in agreement with 
the results found numerically for the linear and zigzag 
configuration in ultracold dipolar gases^. 

Finally, we focus on regime (ii)^ namely the phase 
where the zigzag order found in the framework of the 
classical theory is destroyed as a result of quantum fluc- 
tuations, and the disordered linear chain is found. This 
regime is characterized by a wave function whose width in 
the transverse direction is of the order of the distance be- 
tween the two minima of the potential Vq of Eq. (|^ . see 
for instance the numerical result reported in Ref. Il6l For 
this purpose, we assume that the density along y = a(f> 
is given by ny{y) = 7V|4'+(0)p/a with ^'+(0) given in 
Eq. (gni) and Af = 1/(2 + 2r) giving the right normaliza- 
tion, where F = exp(— (/)q/Z^) is the overlap factor [Eq. 
(1281) ]. The corresponding structure form factor reads 

F + cos(fcj,a</)o) 



Sdis{kx,ky) ~ N^6f,^ 2j^S{kx,ky) 



r + 1 



(56) 

where non-locality is here in the dependence on the pa- 
rameter F, giving the overlap between the two wavefunc- 
tions in the two wells. This result was already obtained 
for the case of a single atom in a double well potential^^, 
and can be understood as an interference pattern aris- 
ing from a photon, which is elastically scattered by one 
of the atoms. In the disordered phase the scattered pho- 
ton will follow two possible pathways, associated with the 
the coherent superposition of the two positions where the 
atom can be found with largest probability. Figure [2Ka) 
displays the function Sdis as a function of ky. The ex- 
istence of sidebands (Fig. [2](a)) is the signature of the 
regime (ii). This demonstrates how one can distinguish 
the various regimes and measure the shift of the critical 
point resulting from quantum tunnelling. 

We also remark that autocorrelation functions of the 
crystal can be measured by means of Ramsey interferom- 
etry, by driving the internal transition of one ion in the 
chain, as proposed in Ref. HI. This method would allow 
one to extract the critical exponents characterizing the 
critical behavior. Methods for revealing quantum tunnel- 
ing between zigzag configurations of few ions have been 
proposed in Ref. 133 . 



VI. CONCLUSIONS 

In this paper, it has been argued that the structural 
phase transition between a chain and a zigzag for ions 
or atoms interacting via a power-law interaction falling 
off as l/r" (where r is the distance and ct > 1) can 
be mapped on an Ising model in a transverse field. In 
this way the classical theory for the phase transition can 
be extended to account for quantum fluctuations. Such 
correspondence was proposed in the past for electronic 
systemgi. Here the mapping is explicitly justified deep in 
the ordered zigzag phase [see Eq. ([50)) ]. Based on sym- 
metry arguments [see Eq. (pS)) ] we postulate that this 
mapping can be extended to the critical regime. The 
number of universality classes in two dimensions is re- 
stricted by conformal field theorj«22. Therefore, we argue 
that identifying the symmetry (Z2 in our case) deter- 
mines the universality class, and implies that it is the 
same as of the one dimensional Ising model in a trans- 
verse field. In other words, we assume that all opera- 
tors obtained in our theory that differ from the ones of 
the conformal field theory of the one dimensional Ising 
model in a transverse field are irrelevant in the sense of 
the renormalization group. To the best of our knowledge, 
there are no existing theoretical tools to verify such state- 
ments. However, this hypothesis is further supported by 
numerical studie o^^i^^ . 

Under these assumptions, we have estimated the quan- 
tum critical point as well as the regime of parameters 
where classically zigzag order would be found, but quan- 
tum tunneling and fiuctuations suppress it and the linear 
chain is recovered. In this calculation it was assumed that 
the effective potential Vq [Eq. obtained in the vicin- 
ity of the critical point according to the classical theory, 
still holds. 

Our analysis allows us to consider possible realizations 
of a suitable system in which the predicted quantum 
phase transition would be observable. In particular for 
trapped ion systems, we identify the physical parameters 
determining the critical point, as well as the experimen- 
tal parameters which are required in order to access this 
regime. Schemes for measuring the transition are pro- 
posed, which are based on photon scattering. 

We remark that this work addresses the two- 
dimensional case, in which the transition to a zigzag may 
only happen in a predetermined plane, and give also the 
full phase diagram including the thermal fluctuations. In 
three dimensions, where a Goldstone mode is predicted 
corresponding to rotations of the plane, in which the 
zigzag is formed around the axis of the chain, one could 
consider an extension of this mapping. On this basis we 
conjecture that the system at the quantum phase tran- 
sition can be mapped to the XY model, and the phase 
transition is of Kosterlitz-Thouless type. 

We flnally note that our theory provides a general 
framework which allows one to study further effects, such 
as dynamics at the quantum critical point. It is worth 
mentioning that here one expects the creation of topolog- 
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ical defects, i.e. domain walls. It would be interesting to 
apply this theory to defect formation, extending the work 
done on the classical system in Ref. [l^, and to explore 
whether these systems allow for the creation of special 
kind of entangled states of the ions. 

To conclude, using a quantum-field theoretical descrip- 
tion we argue that the linear-zigzag instability in two- 
dimensional systems of trapped ions or polar molecules 
can be mapped to the one dimensional Ising model in a 
transverse field. This result demonstrates once more the 
potentialities offered by these systems as quantum simu- 
lators^^i^i^i^, and more generally for quantum techno- 
logical applications. 

Acknowledgments 

We gratefully acknowledge useful discussions with E. 
Altman, E. Demler, J. Eschner, R. Fazio, J. Feinberg, 
M. E. Fisher, Y. Gefen, J. Meyer, M. Raizen and S. 
Sachdev. This work has been partially supported by 
the Israel Science Foundation (ISF), by the US-Israel 
Binational Science Foundation (BSF), by the Minerva 
Center of Nonlinear Physics of Complex Systems, by 
the Shlomo Kaplansky academic chair, by the European 
Commission (IP AQUTE, Strep PICC), by the ESF (EU- 
ROQUAM, CMMC), and by the Spanish Ministerio de 
Ciencia y Innovacion (QOIT, Consolider-Ingenio 2010; 
QNLP, FIS2007-66944; Ramon-y-Cajal). E. S. acknowl- 
edges support from the Ministry of Science and Technol- 
ogy (Grant 3-5792). G. M. acknowledges the German 
Research Council (DEC) for support. 

Appendix A: Semiclassical Energy Eigenvalues for 
the Double-well Potential 

In this appendix the eigenvalues of the Hamiltonian 
Hq [Eq. (HH] are calculated in the semiclassical limit. 



In this limit the term proportional to (j)'^ in the potential 
can be ignored. The energy levels can be found by the 
Bohr-Sommerfeld rule where the action / is quantized: 

/„ = n + const « n (Al) 

for large n. Here, similarly to Eq. (UU, we have used 
units where length scales are normalized by A = y''h/muj 
and energy scales by huj/2. For large energy (E — 
E/{huj/2)) where the semiclassical approximation is rel- 
evant. 



/ = j pd4> = — j \J E - B(l>^d(t) (A2) 

where p is the conjugate momentum of the normalized 
coordinate 0, and B — hg/2'w?uj'^ [see Eq. (03])]. The 
integral is over a closed trajectory of the particle. By 
elementary rescaling one finds: 

I={-in)' C^^fV^'dy. (A3) 

Note that the normalized energies scale as B^^^, as we 
found also for low energies within the approximation 
B >• 1 [see Eq. (Hi])]. The semiclassical energies are 
therefore 

En = B^/^{^y\ (A4) 

The energy spacings increase with energy. In the physical 
units the eigenenergies are given by 
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